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Abstract 

We derive and solve exactly the Dyer-Roeder equation in a Fried- 
man Robertson Walker cosmological model with non zero cosmological 
constant. To take into account non homogeneous distribution of mat- 
ter we use the phenomenological dumpiness parameter a. We propose 
also a general form of an approximate solution which is simple enough 
to be useful in practical applications and sufficiently accurate in the 
interesting range of redshifts. 
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1 Introduction 



Recent determination of distances to high redshift type la supernovae re- 
vealed that the universe is now expanding faster than expected. This ac- 
celerated expansion is linked with a non zero cosmological constant. Since 
gravitational lensing of cosmological objects involves also high redshift galax- 
ies it is interesting to extend the cosmological lens equation to cosmological 
models with non zero cosmological constant. 

Furthermore the mass distribution in the universe at the relevant range of 
redshifts is clearly non homogeneous. In our considerations we will take 
into account the fact that at the scale up to several hundred Mpc matter 
distribution is non homogeneous. Moreover, with the advent of the Next 
Generation Space Telescope (NGT), which will allow imaging of objects at 
z > 10 (( [Barkana fc al.2000| )), it will be possible to check better and better 
the influence of non homogeneities . 

Using the general equations describing propagation of light in a given space- 
time we derive an equation for the angular diameter distance and we extend 
the notion of angular diameter distance to two objects at arbitrary red- 
shifts. The general equation for the angular diameter distance, the so called 
Dyer-Roeder equation, can be exactly solved also in the case when the cos- 
mological constant is different from zero. However, the exact solution is very 
complicated and therefore not useful for practical applications. We propose 
a simple analytic form of an approximate solution. The approximate form 
of the angular diameter distance depends on four arbitrary parameters. We 
fix the values of these parameters by fitting this form to the exact solution. 
The paper is organized as follows: In Section 2 we present the derivation of 
the angular diameter distance in an arbitrary Friedman- Robertson- Walker 
cosmological model. To take into account the non homogeneous distribu- 
tion of matter, following the standard practice, we introduce a constant 
phenomenological parameter d and rewrite the final equation in a few dif- 
ferent forms. In Section 3 we present exact analytic solutions for the angular 
diameter distance, first in a flat universe and then in the general case of ar- 
bitrary spatial curvature. Section 4 is devoted to the approximate solution 
and to discussions of the fitting procedure used to determine the values of 
arbitrary parameters. Discussion of our results is presented in Section 5. 
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2 General considerations 



Let us consider a beam of light emanating from a source S. The hght rays 
propagate along a surface S which is determined by the eikonal equation 

g"^S,,S,^=0 (1) 

A light ray is identified with a null geodesic on S with the tangent vector 
ka = — The light rays in the beam can be described by x°' = x°'{v,y'^) 
where v is an affine parameter and y"" (a=l, 2, 3) are three parameters spec- 
ifying different rays. The tangent vector field to the light ray congruence, 
k" = = — Sja, determines two optical scalars, the expansion 6 and the 
shear a: 

6 = ^A:";, a = k^,pfnf'm'' , (2) 

where m" = "^(^" — ^^") is a complex vector spanning the spacelike 2- 
space (the screen space) orthogonal to k" {k°'fha = 0) (actually the vorticity 
connected with the light beam is zero in all our considerations, therefore in 
our case these two scalars completely characterize the congruence). These 
two optical scalars describe the relative rate of change of an infinitesimal 
area A of the cross section of the beam of light rays and its distortion. In 
particular 

These two optical scalars satisfy the Sachs ( dSachs &: Kristian 19661) ) prop- 
agation equations 

^ + 02 + |a|2 = -li?^^A;"A;^, (4) 

and 

a + 2ea = -^Caf3^sm''k'^m^k\ (5) 

where the dot denotes the derivative with respect to v, Rap is the Ricci 
tensor, and Cap-yS is the Weyl tensor. Equations (^) and (|5|) follow from 
the Ricci identity. We will use these equations to study propagation of light 
in the Friedman-Robertson- Walker (FRW) spacetime. The FRW spacetime 
is conformally flat, so in such spacetimes Ca/3jS = 0. From equation dH) it 
follows that if initially the shear of the null ray congruence is equal to zero 
than it is always zero. Therefore assuming that the light beam emanating 
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from the source S has vanishing shear we can disregard the shear parameter 
altogether. Using (|^ we can rewrite equation (^) in the form 

VA + ^R^pk'^k'^VA = 0. (6) 

An observer moving with the 4-velocity vector u'^ will be associated 
with the light ray circular frequency lo = cu^ka- Different observers will 
assign different frequencies to the same light ray. The shift of frequencies 
as measured by an observer comoving with the sources and an arbitrary 
observer is related to the redshift by 

l + z = — = — n^fca, (7) 

where is the frequency measured by the distant observer. Differentiating 
this equation with respect to the affine parameter v we obtain 

J = -k'^k^u^,^. (8) 




Since the angular diameter distance D is proportional to v A we can 
rewrite equation (P) using D instead of VA and at the same time we replace 
the affine parameter v by the redshift z, we obtain 

"^"^ + ^T^,k-k^D = 0, (9) 

where we used the Einstein equations to replace the Ricci tensor by the 
energy-momentum tensor. 

To relate a solution of @ with the distance it has to satisfy the following 
initial conditions: 

D{z)U=o = 0, 

(10) 

dD{z)^ _ c 
dz ~ Hq 

To be able to use solutions of the Eq. (^) to describe gravitational 
lenses we have to introduce the distance between the source and the lens 
D{zi,Zs), where zi and Zg denote correspondingly the redshift of the lens 
and the source. Let D{zi, Z2) denote the angular diameter distance between 
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a fictitious observer at zi and a source at Z2, of course D{0,z) = D{z). 
Suppose that we know the general solution of equation (P) for D[z) which 
satisfies the initial conditions (0), then the function D{zi,z) defined by 



D{zi,z) 



■il + zi)D{zi)D{z) 



dz' 



D^{z')g{z') 



(11) 



such that 



D{zuz% 



0, 



^D{zi,z)\:,=z^ = sign(z - zi) 



dz 



l + Zi 

9{zi) 



satisfies equation (0), if the function g{z) is a solution of 



d_ 

dz 



\ng{z) 



'dz\ 
dv 



2 ' 



SO 



g{z) = go exp 



zdz . 



(12) 



(13) 



50 is an arbitrary constant of integration, which, without restricting gener- 



ality, we assume to be one. To obtain the equation (12) we inserted (|ll| ) 
into the equation (^) and demanded that it is satisfied. Let us note that the 
Etherington reciprocity relation ( ([Etherington 1933|) ) 



-0(^1, ^^2) _ D{z2,zi) 

l + Zl ~ I + Z2 



(14) 



follows directly from (|Tj 

We are interested in applying equation (^) to find the angular diameter 
distance to objects at high redshifts. Therefore let us consider the standard 
FRW spacetime described by the line element 



ds' 



dt^ - R^{t) 



dr^ 



1 — kr"^ 



(15) 



where k (-1, 0, +1) is the curvature constant, and R{t) is the scale factor. 
If the matter content of the universe can be approximated by dust then, as 
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shown by Friedman, the Einstein field equations for the metric (15) assume 
the form 



^ _ 2 _ ^ Ac^ kc^ 



R _ AttG Ac2 



(16) 



R 

where here the dot denotes derivative with respect to t, H = —, Qm is the 

R 

matter density, and A is the cosmological constant. 

In the FRW spacetime the redshift z is connected with the scale factor 
R 

R{t) by 1 + z = . Differentiating this relation with respect to time we 
R{t) 

obtain 



Let us now return to the description of propagation of a light beam. 
Assuming that the observer is a standard FRW observer, e.g. it is comoving 
with matter, from the equation (^), we obtain 

-^ = il + zfHiz). (17) 
Wo dv 

Introducing a new dimensionless affine parameter w = we trans- 

form equation ([T7| ) into 

Ho— = {l + zfH{z) (18) 

The Friedman equation for H^, Eq. (16), can be rewritten as 

H\z) = Hi (l7„(l + zf + Vtk{l + zf + I^a) (19) 

where -f^o is the present value of the Hubble constant, and 
_ STrGgQ _ kc^ _ Ac2 
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are the present density parameters of matter, curvature and cosmological 
constant respectively, and, in the case p = (no radiation), we have 

nm + ^k + ^A = 1. (20) 

Substituting ([l|) into (||) we get 
dz 



{i + zy^nmii + z)^ + nkii + zy + nA- (21) 

To apply the angular diameter distance to the realistic distribution of 
galaxies it is necessary to take into account local non homogeneities. Unfor- 
tunately so far an acceptable averaging procedure for smoothing out local 
non homogeneities has not been developed (( [Krasinski 1997 )). Therefore fol- 
lowing previous discussions we introduce the phenomenological parameter a 
which describes the influence of local non homogeneities on propagation of 
light ( dPyer fc Roeder 19721) ), (( p:omita et al. 1999| )). With this alteration 
equation (^), in the FRW-dust case, can be rewritten in the form 

dz\^d^D ( d^z\ dD 2, ^ , 

It is customary to measure cosmological distances in units of introducing 
the dimensionless angular diameter distance (the Dyer-Roeder distance) r = 
and using (p^) we finally obtain 



{l + z) + zf + VLk{l + zf + 



dz"^ ^ 



(23) 



7 \ dr 3 

-flmil + zf + 317fc(l + zf + 21^A j ^ + + = 0- 



with the initial conditions: 



r{z)U=o = 0, 

(24) 

dr{z) 



\z=0 



1. 



dz 

This equation can be cast into a different form by using a = as a 
parameter instead of z, we obtain 



(pj. 3 dr 3 

a^i^lm + ^ka + ^lAa^)-r^ - a{-^}m + ^ka)— + -a^mr = 0, (25) 
aa^ 2 da 2 
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or using the cosmic time t as a parameter, the equation ( ^21 ) assumes the 
form 

■^-H{t)— + 47TGag,m{t)r = 0. (26) 

Equation ( p6[) needs some comments: this equation was first introduced 
by Dashevski & Zeldovich (( Dashevski fc Zeldovich 1965| )) (see also Da- 



shevski & Slysh, ( Dashevski &: Slysh 196^ )). More recently Kayser et al. 



(( [Kayser et al.l997 )) have used it to derive an equation similar to our (|25 
In Eq.(|26D the clumpiness parameter a is usually considered as a constant. 
However in the papers by Dashevski & Zeldovich and by Kayser et al. a is 
allowed to vary with time but only the case a = const, is really considered. 
For a discussion of the case in which a depends on z see, for example, the 



paper by Linder (( Linder 198' 



As an example of our procedure let us consider the case when f^A = 
and r^fc = 0, that is a flat Universe so the universe is flat. In this case it is 
easy to see that 

r(z) - + + .27) 

2/3(1 + z)5/4 ' ^^^> 

where /? = — 24q, is the general solution of equation ( p3|) with ap- 

propriate initial conditions. This solution can be found in SEF (( Schneider 



Ehlers &: Falco 1992| )), see their equation (4.56). The g function can be 
easily deduced in the flat A = Friedman-Robertson- Walker cosmological 
model and we obtain 



giz) = {l + zY^nMoZ + l. 

Substituting the expression for g(z) we get the familiar solution D(zi,Zs 
found in SEF. 



3 Exact solutions 



In recent papers Kantowski ( (|Kantowski 1998D ),(( i^antowski fc Kao 2OOO0 ) 
has found the general solution of the Dyer-Roeder equation written in the 
form (p3|). In what follow we use the DR equation in the form (25) and 
obtain the general solution using more direct approach. 
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3.1 The case = 

When Qk = 0, it is possible to divide the Eq.(p5|) by Qm and it becomes : 
where ^ = — — . 

To solve Eq.(|2g) we use the following strategy: First we solve this equa- 
tion in the case when /i = and we look for solutions in the power form a*. 
Inserting this form into (pSl), we obtain: 



2s^ - 5s + 3q = 0, (29) 

which has the solutions : 



s± = t±tV25 - 24a = ^±/3, 
4 4 4 



13 = \\/2b - 24a. 
4 



(30) 



Writing the general solution is terms of z instead of a and imposing the 



initial conditions (11) we recover the solution (|27|). 

In the general case when /i 7^ we look for solutions in the form 
f\ = a^f{x) where x = a^. Inserting this form into ( p8|) after some re- 
arrangements we obtain 

3x{l + ^^^)^+ (2{s + 1)(1 + ^x) - ^) £ + |(s - d)/ = 0. (31) 

This equation can be reduced to the standard hypergeometric equation. The 
general solution of Eq. ([2^) can be written in the form : 



(1 + ^)5/4-' n 1 + 2 



(32) 



' (1 + 2)5/4 Ml 



where Ai,A2 are arbitrary constants and we denoted the solutions by r\ 
to stress that it is the solution of the DR equation in a spacetime with 
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A / 0. Here fs_ and fs^ are solutions of Eq.(^) with s = s+ and s = S- 
correspondingly. The constants Ai and A2 are determined from the initial 
conditions (p!l|), or explicitly: 



dfA, 



dz 



\z=Q 



Aif,Jl)+A2fs_a) = 0, 



-s+Aifs4l) + Ai—fs^\,=o - S-A2fs_il) + 

^23-/5- |2=0 = 1- 

dz 



(33) 



To find the function r\{zi, Z2) we have to solve the equation (|l^), which in 
this case is: 



from which we get: 



9A^ ^ 3 

QA 2a(l + na^) 

1/2 



9A 



[1 + zy + fi 



(34) 
(35) 



Then the general two points {z > zi) solution of Eq.(|28|) is : 
fA{zi,z) = fA(2;i)(l + zi)fA(2;)x 



(36) 



zi I A 



-^{z'){i + z'mi + z'fnm + nAV/^' 



which satisfies the initial conditions: 

rAizi,zi) = 0, 



(37) 



drA x| _ sign(2: - Zl)^/^l^ 
dz ^^i'^^l-=-i - (l + zi)^(l + z)aj7„ + 17A' 



(when Qa = 0, we get exactly the SEF solution (4.53) for = !)• It 
is worth to stress that using properties of the hyper geometric functions 
(/s+)/s_) it turns out that, 

r < fA, (38) 
which is also shown in the following figure. 
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1 2 3 4 5 



Figure 1: The figure shows r{z) for fixed a = 0.9 and for Q\ = 0, Oa = 0.6, 
= 0.65 and Oa = 0.7. These figures have to be read from the bottom 
for increasing values of We see that f < ta as clarified by Eq. (|3^) 

3.2 The case fi^ 7^ 

When the Eq.([25|) can be rewritten in the following form: 



(Pr (5a + I dr 

da? a{a — ai){a — a2){a — 0,2) da 

(39) 

^3^ 1 ^ _ ^ 

2 a^(a — ai)(a — a2)(a — 02) 

where ai,a2,a2 are the roots of the equation 17^, + fi/jO + Q\a'^ = 0, and 
5 = (we are using the symbol r for the Dyer-Roeder distance in this 
more general case) . These roots are : 

V^(-9V^ + V3V4<5=^ + 27^)1/3+ 21/332/3^ 

(1 + iV3)5 



02 



22/33i/3^(_9^+ ^^53 + 27^)1/3 

(1 - i^/3)(-9^ + y3v/4^3 + 27^)^/3 
221/332/3^ 



11 



r 

.5r 




Figure 2: Same as Figl, the z-range is now bigger (0, 100) 



(1 - iV3)S 

~ 22/331/3 ^(_9y^ + ^^/WTrrJiy/^ ~ 

_ (1 + iV3)i-9^ + VS^U^ + 27^)^/3 
221/332/3^ 

where fx = n\/Vtn as before. (02 and are complex conjugate and ai < 0, 
then Eq. ( ^0|) has non singular coefficients for real a > 0). Equation ( |40| ) 
is of Puchsian type (( pnce 1956| ),( [l¥icomi 196l )) with four regular singular 



points plus a regular singular point at infinity. This equation can be put in 
the following form 



(fr I A B C D \ dr 

+ + -1- + 



da^ \ a a — ai a — a2 a — I da 



I fA B C D , 
- - + + H ) r = 0, 

a \ a a — ai a — 02 a — a^, 



(40) 



where the coefficient A, B, C, D, A, C, D are easily found in terms of 
the cosmological parameters ^^a> ^fc and (5, and are given in Appendix 
A. 

The solutions around each of the 4+1 singularity points ( here " 1" 
denotes the singularity point relative to the infinity) are grouped together 
using the so called Riemann-P symbol. 
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ai a2 as oo 
V I 6ii 621 (531 61 a 

^12 ^22 5z2 ^42 ^2 



where the 5ij and 5i are given by 





/f- 


25a 
4 




/f- 


25d 
4 


0, 







-^11 = n-l/-- — (= 2^-), 

2s+), 



522 = 1 + —j^ 77 r, 

ai(a2 - ai)(a2 - as) 

hi = 0, 



'532 = 1 + ^ 



a2(a2 - ai)(a2 - as)' 
hi = 0, 

O42 — i H 7 w T, 

asCas - ai)(a3 - a2) 

<5i = 0, 



<52 = 1, 

and correspond, respectively, to the solutions of the indicial equation relative 
to the finite and infinite regular singular points of the equation. We do 
not discuss further properties of these solutions because unfortunately they 
cannot be given in an explicit analytical form. In the next section we will 
give an approximate analytical expression for the solution of equation (|25|) . 
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4 The approximate lens equation 



The exact solutions presented in the previous section are comphcated and 
difficult to use in practical applications. Therefore we want to find an an- 
alytical approximate expression for r{z,il.jn,Q\,Qk,0() and for the function 
X appearing in the cosmological lens equation. Following SEF let us briefly 
recall the basic equations used in the derivation of the cosmological lens 
equation and introduce appropriate notation. The time delay between dif- 
ferent light rays reaching the observer is: 



cAt = (1 + Zd) 



DdDs 
Dds 



13) - V'(C) \ + const 



(41) 



where D^, Dg and D^g are correspondingly the angular diameter distances 
to deflector, source and the angular diameter distance between deflector and 
source. 

The first term in the bracket represents the geometrical time delay and 
the second one is connected with the non homogeneous distribution of mat- 
ter. To describe this effect we use a perturbed metric in the form : 



where r] is the conformal time, U is the gravitational potential of the deflec- 
tor, and da"^ is the spatial line element. The first term in Eq.(^l|) can be 
rewritten as : 

DdDs 



(1 + Zd) 



D. 



ds 



ixizd) - x{z.s)y 



(42) 



where the function x is given by : 



X{z,flm,^A,^k,a) 



dz 



'9{z) 



exp 



d'^z 



dz\' 
dX J 



r-^dz, (43) 



so it is connected with the general solution of Eq.(23) and the general ex- 
pression for g{z) (see Eq.|25D. Since the metric is conformally stationary, 
from the Fermat principle, we get the cosmological lens equation : 



c-tio 86 



(44) 
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where Rg is the Schwarzschild radius of the deflector. Denoting = 0^9, 
r] = Usp, a(4j = -j^ -, we obtain 

r/=^e-I?ds5(e) , (45) 

which is formally identical to the lens equation obtained for z << 1. 

Let us first discuss the interesting range of the parameter space. From 
the available observations and theoretical considerations it follows that the 
radius of curvature of the universe is very large or equivalently that dimen- 
sionless curvature parameter 0,^ is very small. In our considerations we 
assume that Ififcl < 0.05. Unfortunately no reliable estimates of a exist, but 
in the considered range of redshifts ze(0, 100), we consider the range of a 
£ [0.3,1] (not completely clumpy Universe). We also allow different values 
of JIa and with special attention we treat the case f^A = 0.65 which seems 
to be singled out by o bservations ((p<Cochanek 19961), ( [Kochanek 1996bD , 
(iKochanek fc al.l998| ), (|Perlmutter 19971) , ([Tomita 1999|) ). The density pa- 
rameters satisfy the constraint + + ^fc = 1- 

In the considered range of redshifts and of the other cosmological pa- 
rameters we propose to use the following form of the approximate function 
for r 

r(z) = ^ (46) 

where the parameters di , ^2 , da , (i4 depend on the cosmological parameters 
and they are determined by fitting Eq.(^6|) to the corresponding numeri- 
cal solution. To optimize the fit we use a non linear regression procedure 
available in Mathematica 4-0- 

We use the same procedure to find the approximate function x{z) which 
we take in the form 

(ei + e2Z + esz^) 

where the parameters ei, 62, 63 are determined by fitting xi^) to numerically 
obtained exact values. 

We are interested in finding an analytical approximation of the Dyer- 
Roeder distances for cosmological values of redshifts (say, for example, z > 
0.05); in this region the polynomial, which appears in the denominator of 
the approximate function xi^) does not have zeros. In Fig. 5 - 6 we show 
the exact angular diameter distance and the fitted approximate relation for 
the same value of a and other cosmological parameters. (The calculations 
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r 

0.5r 




z 

20 40 60 80 100 

Figure 3: r{z) for different values of a, and for fixed CIa = 0.65 is shown. 
The values of a are a = 0.9, a = 0.8, a = 0.7 respectively, a decreases from 
the bottom curve up. 



have been done for fixed a, = 0.65 and Qk = 0.) In Fig.7 we see that 
the approximate angular diameter distance reproduces the exact curve with 
a very good accuracy, the error is not larger than ~ 1%, for z in the range 
(0, 50) and it is still good for larger z in the range (50, 100). In Fig. 8 we 
show exact numerically obtained function x(z) and the fitted approximate 
relation for the same value of a and other cosmological parameters. 

Finally we have analyzed how the typical maximum present in the angu- 
lar diameter distance depends on the two cosmological parameters (i^Ai ck)- 
Actually we perform this analysis both for the maximum of that distance 
rmax and for the z- value for which that maximum occurs Zmax- For Qk = 
and J7a = const., Vmax depends on a as: 

rmax{a, Oa = const., = 0) = 

(48) 

= exp[Co + Ci(5^a) + C2(5^a)' + CsI^^a)'] 
If we fix a = const, and consider the dependence of rmax on Oa, we find: 
rmax{a = const., Oa, Ofe = 0) = 
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0.5 
0.4 
0.3 
0.2 
0.1 



1 2 3 4 5 

Figure 4: Same as in Fig. 3; the z- range is now(0,5) 

(49) 

= exp[eo + eiinA) + e2(5^A)^ + ^sI^^a)^], 

where the coefficients Si depend on the dumpiness parameter a and on O^. 
We can also study how Zmax depends on a and J7a- For a = const., = 0, 
we have that 

Zmaxi^ = const., QA,^k = 0) = 

(50) 

= exp[ro + ri(!^A) + T2{nAf + Tsin^f + u{nA)^ + M^Af]. 

The coefficients Tj depend on the dumpiness parameter a and on fi^. For 
r^A = const., Qfc = 0, we find 

Zmaxia,^A = const. ,Qk = 0) = 

(51) 

= exp[7o + 7id + 72(0;)^ + 73(0)^ + 74(0)^ + 75(0;)^]. 

where the 7^ coefficients depend on and the parameters I^a and . From Eq. 
( |5l| ) it follows that for Oa = 0.65 and a = 0.9, we get Zmax = 1-62. When 
Oa = 0, Ofc = 0, as it follows from Eq. (p7|) Zmax = 1-25. It is interesting 
to note that Zmaxi^A = 0, fi^ = 0, d = 0.9) < Zmaxi^A = 0.65, fifc = 0, d = 
0.9), which shows how relevant is the role of the cosmological constant. 
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Figure 5: The numerical solution of the Dyer-Roeder equation and our ap- 
proximated function for fixed I^a = 0.7 and a = 0.9 are shown. As it is 
apparent, the approximating function works very well. 

5 Conclusion 



In this paper we have considered the cosmological lens equation in the case 
when the cosmological constant A 7^ and the curvature of space is different 
for zero (A: 7^ 0). We have included the effects of non homogeneous distri- 
bution of matter which are described by a phenomenological parameter a. 
Unfortunately at the moment there are no generally accepted models that 
describe the distribution of baryonic and dark matter at high redshifts and 
therefore the influence of non homogeneities of matter distribution can be 
included only at this approximate level. Following the standard procedure 
(( [Schneider, Ehlers fc Falco 1992|) ) we use the Dyer-Roeder distance to find 
the distance between two objects with redshifts zi and Z2- 

To find the general solution we slightly enlarged the parameter space by 
considering the non zero cosmological constant, so the cosmological model 
is described by the following parameters Q\ ^ 0, / 0, d G [0.3, 1] (not 
completely clumpy Universe), p = (no radiation). The solution that we 
have found is a functional combination of the solution given in SEF for the 
case when Q\ = 0, $7^ = and two hypergeometric functions. Fig. 1 shows 
that in the flat universe the Dyer-Roeder distance increases for increasing 
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Figure 6: Same as in Fig. 5, the z-range is now (0, 5.) 



values of the cosmological constant. In the next step we allow the curvature 

of space to be different from zero. This further complicates the Dyer-Roeder 
equation which becomes of a Fuchsian type with 4 regular singular points 
and one regular singular point at infinity. The general form of this kind of 
ordinary differential equations is given in terms of the P-Riemann symbol. 

The obtained exact solution is so complicated that it is practically use- 
less in practical applications. Therefore we have looked for an approximate 
analytic solution simple enough to be used in many applications and at the 
same time sufficiently accurate at least in the interesting range of redshifts. 
The proposed form of the approximate solution of the Dyer-Roeder equation 
depends on four arbitrary parameters. To fix values of these parameters we 
fit the approximate solution to the exact one with the help of a non linear 
regression method (see Mathematica 4-0) ■ Of course, the parameters de- 
pend on the density parameters, but their values for different cosmological 
parameters can be tabulated. Following SEF we have also found the func- 
tion X which appears in the expression for time delays as well as in the lens 
equation. We have also proposed approximate analytical form of the func- 
tion X which contains three arbitrary parameters. To find values of these 
parameters we use the same method as above. 

We have also studied how the redshift corresponding to the maximal 
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angular diameter distance depends on the basic parameters determining the 
cosmological model. 

Let us now consider the following three combinations of the DR distance: 

Dls ^ Hq ris 
Dos c ros' 

DolDlS _ Hp roLTLS 

Dos c ros 

DolDos Hq roiros 
o = . 

Dls c rLS 

We have selected these combinations because of the role they play, respec- 
tively, in bending of light, lensing statistics, and time delay. Asada ((Asada 
19971 ), ( lAsada 199^ ), ( |Asada 1998|) ) have found that: 



-^^("i) < -^^("2), for ai < 02, 
Dos Dos 

DolDls , N ^ DolDls , N c 

— f; (ai)< — 7^ (a2), for qi < a2, 

Dos Dos 

DolDos , s ^ DolDos , X r 

— (ai) > — (02), for ai < a2- 

Dls Dls 
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Figure 8: This figure shows how our approximate analytical form for x 
interpolates the numerical values. The values of non normalized chi-squared 
of this regression is = 0.013. The independent variable is z-lO. 



Using our approximate relations for r and x we obtain that: 
Dls 



Dos 



DqlDls 
Dos 



DolDos 
Dls 



(1 + ZL)r{zL){x{zL) - X{zs)) 



^di + d2{d2. + di/zL + zlY) (ei + e2^L + e^z]) 



1 



(62 + 62^5 + es^s) 



{I + ZL)r^{zL){x{zL) - X{zs)) 



(1 + zl) 



1 



{di + d2{dz + d^/zL + zlY) (ei + e2ZL + e^zD 
1 



(ei + 62^5 + es^l) 
1 

{I + zl){x{zl) - X{zs)) 
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Using the fitted values for the parameters (di, (i2, da, ^4) and {ei, 62,63) we 
confirm results of Asada for a in the considered range 1 > d > 0.3. 
Finally we like to stress that from our analysis it follows that variations 
in the angular diameter distance caused by the presence of cosmological 
constant are quite similar to variations due to changes in the value of d. 




Figure 9: r(z) for = and d = 0.6 (upper curve) and JIa = 0.7 and 
d = 1 (bottom curve) . This figure shows that non homogeneous distribution 
of matter can mimic the effect of non zero cosmological constant. 

Actually, in Fig.9, we consider two plots corresponding to the values 
Q\ = and d / and JIa 7^ and d = 1 and we see that the d parameter 
could mimic the effect of a non zero cosmological constant. This is an impor- 
tant observation in view of the recent observational results concerning the 
non zero value of the cosmological constant (( [Kochanek 199^ ), ([Kochanek 
&: al. 19981 ), ( [Kochanek 1996b| ), ( PerlmutterT997D ). In the fut ure we would 
like to extend our work and include also the radiation density parameter. 
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Appendix A 



2aia2a3 ' 

^ 3a^m ^ —cx^ic^s — 0.4) + cti{ai — Q4) — Q8x(ai — 03) 

2ag {-ae{a3 - 0:4) + 0:7(02 - 04) - 08(0:2 - "3)) ' 

C = -{A + B + D), 

^ _ {A{ai - 03) + ■B(oi2 - as)) 
(03 - 04) 

The coefficients A, B, C, D are 



2 010203 '' 



2o;9(a30!6 — 040207 + o;20;8 — Q:2Q!8) 
^ 3aiO!7^^^ - 30:407^^^ + Saias^m _^ 

(209(0:306 — 0;40!2Q!7 + 0:208 — O2O8)) 

_l 2o709^^fc - 2a8Q:9^fc 

(209(0306 — 040207 + 02O8 — 02O8)) ' 

~ —0205^2^ + 80405$! 

m 

209(0306 — 0405 — O2O7 + 0407 + O2O8 — O3O8) 
SoiOe^^ + 30206$!m 



(209(0306 — 0406 — 0207 + 0407 + 0208 — O3O8)) 
SoiOsJ^m - 302a8$!m " 2a6a90jfc + 20809^^ 



+ 



(209(0306 — 0406 — 0207 + 0407 + 0208 — 0308)) ' 

^ 30205^m - 30305 -Oyn + 30306^^ ^ 

2o9(— O3O6 + O4O6 + 0207 — O4O7 — 0208 + O3O8) 

^ 30i07Qm - 302a7^m - 20608^m " 2a6Q:9^fc 
(2o9(— O3O6 + 0405 + O2O7 — 0407 — O2O8 + O3O8)) 
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{2ag{—asa6 + 040^ + 020:7 — 0:40:7 — 02O8 + 0303)) 
The functions Oj being: 

oi = aia2 + aias + 0203, 

(X2 = (12(13, 
"3 = 0103, 
Q4 = 0,1(12, 
05 = ai + 02 + 03) 

ae = 02 + 03, 
07 = ai + 03, 
as = ai + 02, 
og = 010203. 
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7 Appendix B 



We give the form of the coefficient appearing in equations (|4^) and ( p7| ) , 
in function of the cosmological parameter. We show the dependence of the 
coefficient on the dumpiness parameters a, fixing Q\ = 0.65 

di = exp[adi + bdia + Cdio^ + ddidr"], 
d2 = exp [ad2 + bd2a + Cd2a'^ + dd2a^], 
ds = exp [ads + 6^3" + Cds"^ + c^ds"^], 

(i4 = exp [ad4 + bdidi + CdA^c? + ddA,c?\. 
Dependence of di functions from JIa with q = 0.8 

dx = exp \pd\ + qdx^K + ?'di(^^A)^ + Sdi(^^A)^], 

d2 = exp (prf2 + qd2^K + rd2{^kf + Sd2{^kf), 

(is = exp {pdz + Qds^A + Tdsi^Af + Sdsi^hf), 

di = exp IpdA + QdA^A + rd4(^^A)^ + Sd4(^^A)^], 
being Pn{^}\) a polynomial of n-degree in ^}\. 

We see below the dependence of the functions from a with = 0.65 

ei = TTiel + nel O + Ogl + PelCe^; 

62 = me2 + ne2a + Oe2a^ +Pe20?, 

63 = + riesa. 

Here, we show the dependence of the functions from OA,with a = 0.9, 
and r^fc = 0: 
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ei = m~ei + n~ei U\ + 0~el ^\ + Pel^^A, 

62 = rne2 + n~e2 + Oe2 ^\ + Pe2 

63 = Pg(riA)/(exp [mesf^A + Wes]), 

being Pq{^K) a polynomial of q-degree in f^A depending on the value of k 
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